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Abstract 

The dipping X-ray source 4U 1822—37 was observed by Suzaku on 2006 Octrober 
20 for a net exposure of 37 ks. The source was detected with the XIS at a 1-10 keV 
flux of 5.5 xlO” 10 erg cm -2 s -1 , and with the HXD (HXD-PIN) at a 10-50 keV flux of 
8.9xlO“ 10 erg cm -2 s _1 . With HXD-PIN, the pulsation was detected at a barycentric 
period of 0.592437 s, and its change rate was reconfirmed as — 2.43 xlO -12 s s _1 . The 
1-50 keV spectra of 4U 1822-37 were found to be very similar to those of Her X- 
1 in the slopes, cutoff and iron lines. Three iron lines (Fe Ka, Fe XXV, and Fe 
XXVI) were detected, on top of a 1-50 keV continuum that is described by an NPEX 
model plus a soft blackbody. In addition, a cyclotron resonance scattering feature 
was detected significantly (>99% confidence), at an energy of 33±2 keV with a depth 
of OHloJ. Therefore, the neutron star in this source is concluded to have a strong 
magnetic field of 2.8 xlO 12 G. Further assuming that the source has a relatively high 
intrinsic luminosity of several times 10 3 ' erg s _1 , its spectral and timing properties 
are consistently explained. 
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1. Introduction 


Most of low- mass X-ray binaries (LMXBs), namely X-ray emitting binaries with Roche- 
lobe filling low mass stars, are considered to involve neutron stars (NSs) with low magnetic field 
strengths, B <10 10 G. They make a contrast to high-mass X-ray binaries (HMXBs), mostly 
containing NSs with high magnetic field strengths (B ~10 12 G) which capture stellar winds 
from their companions. In fact, among ~50 known LMXBs, only three are known to have NSs 
with strong magnetic fields; Her X-l, 4U 1626—67 and GX 1+4 (Nagase 1989). These distinct 
combinations of the mass-donating and mass accreting components may be generally interpreted 
as population effects, that older NSs have weaker fields. However, it is not necessarily clear 
whether NSs gradually lose their magnetic fields (e.g., Chanmugam 1992). Then, it is worth 
while searching for other LMXBs that involve NSs with strong magnetic fields. 

Evidently, the magnetic field affects accretion mechanisms in binaries. When the NS has 
strong magnetic fields, the accreting matter is funneled onto the two magnetic poles, leading to 
strong pulsations, and the production of very hard X-ray spectra which is often accompanied 
by cyclotron resonance scattering features (CRSFs; Makishima et al. 1999). If, in contrast, the 
NS is weakly magnetized, an accretion disk is considered to extend down to vicinity of the NS, 
and the emergent X-ray spectra will exhibit characteristic bimodel behavior between so-called 
soft state and hard state (e.g. Lin et al. 2007; Sakurai et al. 2012). These spectral properties, 
together with the presence/absence of X-ray pulsations, will conversely allow us to tell whether 
the NS in a mass exchanging binary is strongly or weakly magnetized. 

In an attempt of looking for LMXBs that involve strongly magnetized NSs, the present 
paper focuses on the dipping X-ray binary 4U 1822-37, located at an estimated distance of 2.5 
kpc (Mason et al. 1982; Cowley et al. 2003). The companion star in 4U 1822—37 is estimated to 
have a mass of 0.44 - 0.56 M 0 (Munoz-Darias et al. 2005), with M 0 being the Solar mass, and 
hence this system is clarified as an LMXB. X-ray and optical light curves of this source both 
show intensity modulations and dips synchronize with its orbital period, P OT b ~ 5.7 hr (Mason 
et al. 1982; Burderi et al. 2010). These effects are attributed to either gradual occultation (by 
the companion star) of a large X-ray scattering corona above the accretion disk (e.g., Iaria et al. 
2001), or variable attenuation of direct X-rays from the NS by ionized humps on the disk (e.g., 
Diaz Trigo et al. 2006). From dip properties, the orbital inclination is constrained between 76° 
and 84° (Hellicr & Mason 1989; Heinz & Nowak 2001). 

Although 4U 1822—37 is classified as an LMXB, its X-ray properties appear rather 
different from those of other dipping LMXBs. First, it shows clear X-ray pulsations with 
a period of P s = 0.5924 s which is much longer than those of other typical LMXBs, a few 
milliseconds (Wijnands 2004) if ever detected. Figure 1 is so-called Corbet diagram (Corbet 
1984), which displays neutron-star binaries on a plane of orbital period and rotation period. Like 
4U 1626—67 and Her X-l, 4U 1822—37 is thus located on this plot between typical LMXBs and 
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Fig. 1. The Corbet diagram for neutron-star binaries of known rotation periods. Circles are LMXBs and 
triangle HMXBs. 

HMXBs. Combining the estimated X-ray luminosity (10 36 ~ 38 erg s _1 ) with P s and its change 
rate, P s = —2.85 x 10” 12 s s _1 , the magnetic held strength of the NS in this source was actually 
estimated as B = 10 8 ~ 16 G (Jonker & van der Klis 2001). Although the error is very large, it 
is noteworthy that rather high held strength are allowed. 

Second, this object exhibits much harder spectra together with a lower cutoff energy 
than typical LMXBs in the hard state. From the measured cutoff energy (~ 6 keV), Parmar et 
al. (2000) obtained an estimate as B ~ 10 12 G. Third, Torrejon et al. (2010), using the Chandra 
HETGS, detected from this source narrow iron K a and Kg lines, which are generally absent or 
much weaker in other LMXBs (Cackett et al. 2008). 

In the present study, we examine whether the NS in 4U 1822—37 has B 10 12 G or 
not. For this purpose, it is necessary to obtain spectra with a good energy resolution and a 
wide energy coverage. The fifth Japanese X-ray satellite Suzaku (Mitsuda et al. 2007), carrying 
onboard the XIS (Koyama et al. 2007) and the HXD (Takahashi et al. 2007), is best suited for 
the requirements. We hence utilized an archival Suzaku data set of 4U 1822—37, and study 
its spectral and pulsation properties. The mechanism of the periodic X-ray and optical dips, 
though interesting, is beyond the scope of the present paper. 

2. Observations and Data Reduction 

With Suzaku, 4U 1822—37 was observed for an exposure 37 ks on 2006 October 20 (ID 
401051010). In this observation, the target was placed at the XIS nominal position. Of the four 
XIS cameras, XIS0 was operated in the 1/4 window mode, and the others in the full window 
mode. The HXD was operated in the normal mode. 

Using HEAsoft ver 6.12, we analyzed the data from XIS0, XIS2, XIS3 and HXD-PIN. In 
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this study, we did not use XIS1 which is the back-illuminated camera, because it has a higher 
background and is subject to larger calibration uncertainties than the other XIS cameras. As 
the average counts rate of the XIS, ~ 10 cts s -1 , was high enough to cause pile up under the 
full window mode, the XIS events were accumulated over an annular region with the inner and 
outer radii of O'. 5 and 2'.0, respectively. This allowed us to reduce pile up effects to within 1% 
(Yamada et al. 2012). The response and arf hies of the XIS were generated using xisrmfgen 
and xissiraarfgen, respectively. The HXD-PIN data were analyzed using a “tuned” non X-ray 
background hie and the epoch2 response hie, both officially released by the HXD team. We 
did not use the HXD-GSO data because the source was undetectable therein. 

The present Suzaku data were already analyzed by Iaria et al. (2011). Combined with 
previous data, they determined the orbital period to be P or b = 20054.2049 s = 5.57 hr, but did 
not report on pulse detection, or spectral analysis. 


3. Results 


3.1. Timing analysis 

Figure 2 shows background-subtracted light curves of 4U 1822—37 obtained with the 
XIS (1 to 10 keV) and HXD-PIN (15 to 60 keV), together with the HXD-PIN vs. XIS hardness 
ratio. The gross exposure of this observation (~90 ks) covered about four orbital cycles, in 
which we detected about 4 dips in both the XIS and HXD-PIN bands. However, unlike the 
cases of many other dipping sources (e.g. XB 1916—053, XB 1323—619 and EXO 0748—176), 
no X-ray bursts were detected. The same light curves, folded at P or b=20054.2049 s, are shown 
in figure 3. They are consistent with those of previous studies (e.g., Iaria et al. 2008; Iaria et 
al. 2011). From the bottom panel of figure 3, the hardness ratio is observed to decrease during 
the dips. 

We also tried to detect the pulsation at a period ~0.59 s (Jonker & van der Klis 2001), 
using only the HXD-PIN data because of the low time resolution of the XIS. After applying 
barycentric corrections to the individual HXD-PIN events, we further corrected the event arrive 
times for the expected orbital delay, At, in 4U 1822—37. This At is calculated with the pulsar’s 
semi-major axis a and the inclination i as 
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where <% (0 A 4>o A 1) is the initial orbital phase. We fixed asini at 1.006 lt-s, after an 
accurate measurement by Jonker & van der Klis (2001), and chose 0 O — 0.445 so that At 
becomes maximum at the observed X-ray dips. After these corrections, we calculated a 15-40 
keV periodogram. As shown in figure 4(a), the pulsation was detected with a high (> 99% 
confidence) significance at a period of P s = 0.5924337 ± 0.000001 s. The pulse profile folded 
at P s is presented in figure 4(b). As already reported (Jonker & van der Klis 2001), the pulse 
fraction is rather small, ~ ±5% in the relative peak amplitude. 
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Fig. 2. Background subtracted light curves of 4U 1822—37 with 128 s bins, obtained with XIS0+XIS3 
(top panel; 1-10 keV) and HXD-PIN (middle panel; 15-60 keV). The bottom panel shows the HXD-PIN 
vs XIS hardness ratio. 



Fig. 3. The same light curves of 4U 1822—37 as in Hgure2, folded at the orbital period of 20054 s. The 
three panels correspond to those of figure 2 . 

In figure 5, the measured value of P s is compared with previous pulse-period measure- 
ments. Over the past ~ 6 years, the object has thus been spinning up monotonically with an 
approximately constant rate of P s = — 2.43±0.05xl0” 12 s s -1 , or P s / P s = 6.7 kyr. This value of 
P s reconfirms the previous measurements by Jain et al. (2010). 
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Fig. 4. (a) : The periodogram of 4U 1822—37 with the 15-40 keV HXD-PIN data. Abscissa shows period 
difference from a fiducial value of 0.5924337 s. (b) The 15-40 keV pulse profiles, folded at the fiducial 
period in (a). The background is included, at a level of 15% of the total counts. 



Fig. 5. The long-term history of P s of 4U 1822—37 based on Jain et al. (2010) and the present work. The 
triangle indicates the present result. 
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3.2. Spectral analysis 

Figure 6(a) shows 1-50 keV spectra of 4U 1822—37 obtained with the XIS and HXD- 
PIN. The same spectra, normalized to a power-law model of photon index 2.0, are compared 
in figure 6(b) with those of another dipping source EXO 0748—676 and the X-ray pulsar Her 
X-l. The spectral shape of 4U 1822—37 is similar to that of Her X— 1 while different from that 
of EXO 0748—676, in the following two points. First, the spectrum of 4U 1822—37 shows a 
very hard slope (with a power-law photon index of T ~1) below 10 keV and a steep cutoff at 
~ 15 keV, in resemblance to the Her X-l spectrum, while that of EXO 0748—676 exhibits a 
softer (T ~2) slope without clear high-energy cutoff. Second, complex Fe K -a lines are seen at 
~ 6.4 keV in the spectra of 4U 1822—37 and Her X-l, while they are absent in the case of EXO 
0748—676. These suggest that the NS in 4U 1822—37 has strong magnetic fields, like that in 
Her X— 1 (Enoto et al. 2008). 

To quantify the XIS and HXD-PIN spectra of 4U 1822—37, we fitted them jointly with 
an absorbed cutoff-power-law (CutoffPL) model, including three narrow Gaussians to represent 
Fe-Ka, Fe XXV Ka, and Fe XXVI K a lines, at ~6.4, ~6.7 and ~6.9 keV, respectively. The 
photon index T, the cutoff energy E cut , and normalization of CutoffPL were left free, as well 
as the column density of absorption. Similarly, the three Gaussians were allowed to have free 
center energies, free normalizations, and free widths. The HXD vs. XIS cross normalization 
was fixed at 1.16 (Kokubun et al. 2007). However, a large data excess below 3 keV made 
the fit unacceptable, with a reduced chi-square of xl >3.0 for u=2S7 degrees of freedom. To 
account for this “soft excess” feature, we chose a blackbody (BB), and fitted the data with 
wabsx (BB+CutoffPL+3gaussians) . A BB with a temperature of ~ 0.1 keV improved the fit 
to x^=l-21 (zz= 284). However, as shown in figure 7(b), there still remained negative resid- 
uals at ~ 35 keV. Since they are suggestive of a CRSF, we employed a model of the form 
of wabsx (BB+CutoffPL+3gaussians) xcyclabs, where cyclabs represents the cyclotron res- 
onance absorption factor given, e.g., in Makishima et al. (1999). By allowing the resonance 
energy, width, and the depth of cyclabs to vary freely, the fit has been improved to xt = 1-14 
for i^=280. As summarized in table 1, the center energy of CRSF was obtained as -Ecrsf=33±2 
keV. 

To make the model more physical, we once excluded the cyclabs factor, and tried a 
thermal Comptonizetion (CompTT; e.g., Titarchuk & Lyubarskij 1995) and a “negative-positive 
power- law with exponential cutoff’ (NPEX; e.g., Makishima et al. 1999) continua instead of 
CutoffPL, while retaining the BB component and the three Gaussians. As summarized in 
table 1, the CompTT continuum was less successful than that with CutoffPL, while the NPEX 
continuum was better. As shown in figure 7(e), however, the BB+NPEX continuum again left 
the negative residuals at ~35 keV like the case of CutoffPL. We hence revived the cyclabs 
factor, and fitted the data with wabsx (BB+NPEX+3gaussian) xcyclabs. The fit has become 
fully acceptable, x^=1.13 with z/=279. 
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Fig. 6. (a) Background subtracted spectra of 4U 1822—37. The inset shows an expanded view around the 
iron lines, (b) Suzaku spectra of 4U 1822—37 (black; same as panel a), Her X-l (blue), and EXO0748— 676 
(red), all normalized by a common power-law of photon index 2.0. 

The obtained cyclabs parameters, shown in table 1, are not much different from those 
with the CutoffPL continuum case, except that the width has now been constrained as 5.0l|;g. 
While the error ranges in table 1 refer to 90% confidence limit, the cyclabs depth still remains 
positive, D >0.24, if we employ more conservatively 99% limit (i.e. , Ay 2 =6.63 for a single 
parameter). Since the inclusion of the cyclabs factor improved the fit by Ay 2 =— 13.7 for 
Az/=— 3, or Ay 2 /Azz = 4.56, an F-test indicates that the fit improvement is significant at a 
confidence level of 99%. Using the gabs model for the CRSF instead of the cyclabs model 
gave nearly the same xl and a consistent value of Ecrsf From these results, we can claim the 
detection of a CRSF. 

The 1-50 keV flux of this source becomes 1.5 xlO -9 erg s' 1 cm -2 , and the luminosity is 
Lx =1.0xl0 36 erg s -1 at the distance of 2.5 kpc (Mason et al. 1982). These are consistent with 
previous reports (Mason et al. 1982; Iaria et al. 2001). 

4. Discussion 

4-1. Spectral shape 

In figure 6(b), the spectra of 4U 1822—37 were visually compared with those of Her 
X-l and EXO 0748—676. For doing this more quantitatively, we fitted the CutoffPL model 
to Suzaku spectra of other hard-state LMXBs (Aql X— 1, EXO 0748—676, XB 1323—619, 4U 
1636+536 and 4U 1608—52), and accreting pulsars which have strong magnetic fields (4U 
1626-67, Cen X-3, SMC X-l, Her X-l, GRO J1008-57, 1A 1118-61, 1A 0535+26 and LMC 
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Fig. 7. (a) vFu spectra of 4U 1822—37, fitted with wabsx (BB+NPEX+3x Gaussian) xcyclabs. 
The softer and harder continuum components represent the negative and positive 
power-laws of the NPEX model, with a common exponential cutoff factor. (b) 

Fit residuals with the wabsx (BB+CutofFPL+3x Gaussian) model. (c) Those with 

wabsx (BB+CutoffPL+3x Gaussian) xcyclabs. (d) Those with wabsx (BB+CompTT+3xGaussian). (e) 
Those with wabsx (BB+NPEX+3x Gaussian), (f) Those with wabsx (BB+NPEX+3x Gaussian) xcyclabs. 
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Fig. 8. Photon indices and cutoff energies of the Suzaku spectra of representative neutron-star binaries. 
Circles represent hard state LMXBs, squares accreting pulsars, and the triangle 4U 1822—37. 
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X-4). Figure 8 summarized those fit results, on a plane of T vs E cut . The hard state LMXBs 
typically have softer (T ~1.7) slopes and higher cutoff energies (-Ecut > 50 keV), while the 
accreting pulsars show harder (T ~0.6) slopes together with lower cutoff energies ( E ~10 keV). 
Clearly, 4U 1822—37 is more similar to the accreting pulsars than to the hard-state LMXBs. 
In terms of the continuum slope, 4U 1822—37 is thus inferred to have common characteristics 
of accreting NSs with B ~ 10 12 G. 

4-2. The Cyclotron Resonance Scattering Feature 

As argued in Makishima et al. (1999), the high-energy cutoff of the spectra of accreting 
pulsars, which are often steeper than a thermal rollover as pointed out by Tanaka (1986), 
is presumably caused by the presence of a CRSF, which is empirically thought to appear at 
energies of Ecrsf ~ 1-5 x E^p t . Then the NPEX fit results, -E cut =4.8-6.2 keV, predict a CRSF 
to appear at -E C r SF =16-23 keV. This argument was already employed by Jonker & van der Klis 
(2001) to suggest B ~ (1 — 5) x 10 12 G. 

Indeed, as expected, we have detected a CRSF at E GR SF ~ 33 keV with a high signifi- 
cance. Using the basic relation of B=(l + z ) x (-Ecrsf/H-6) x 10 12 G (e.g., Makishima et al. 
1999), where £ is the gravitational redshift of the NS, we obtain 5=(2.8±0.2) xlO 12 G assum- 
ing z= 0 for simplicity. This gives the most convincing evidence that the NS in 4U 1822—37 is 
strongly magnetized. 

Just to make the CRSF detection more convincing, let us examine the derived CRSF 
parameters. As already confirmed, the resonance energy -Ecrsf is consistent with the continuum 
shape. The depth, D ~0.4, is reasonable, in comparison with typical values of 11=0.1-1.7 found 
in other objects. In addition, IU=2-10 keV=(0.06-0.30)E'cRSF agrees, within errors, with the 
general scaling of IU=(0.27-0.50)-Ecrsf found by Makishima et al. (1999). Thus, the present 
CRSF interpretation of the ~33 keV spectral feature is considered reasonable. 

4-3. Spin period and Spin-up rate 

The present HXD data yielded P s = 0.5924337 s and P s = -2.43 xl0~ 12 s s -1 . When 
these values, together with 5=2.8 xlO 12 G, are substituted into the accretion torque formula 
by Ghosh & Lamb (1978), namely 

-P = 1.84 x 10" 12 x B 12 x ( PL% 7 ) 2 s ■ s" 1 (2) 

where L37 is the luminosity in units of 10 37 erg s” 1 and B 12 is the magnetic field in 10 12 G, we 
obtain L~ 3 x 10 37 erg s 1 . 

This luminosity is much higher than the value of ~1.0xl0 36 erg s^ 1 derived from the 
observed flux and an assumed distance of 2.5 kpc (Mason et al. 1982). One possible cause 
of this discrepancy is that the object is in really located at ~ 10 kpc distance, instead of the 
2.5 kpc which is based on some assumptions. An alternative possibility, already pointed out 
previously (White & Holt 1982; Parmar et al. 2000), is that the object appears unusually X-ray 
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faint due, e.g., to X-ray obscuration by some ionized materials on the accretion disk. Indeed, 
its X-ray to optical luminosity ratio of ~20 for Lx ~lxl0 36 erg s" 1 (Mason et al. 1982) is much 
lower than a typical value of ~500 for LMXBs (van Paradijs & McClintock 1995). The ratio 
of 4U 1822—37 will increase to ~600 if we employ Lx ~3xl0 37 erg s 1 . 

4-4- I ron lines 

Another interesting discussion may be performed on the iron lines in the spectra. 
Generally, LMXBs have weak, sometimes broad (Cackett et al. 2008) iron lines with small 
equivalent width (EWs); e.g., 20 ± 15 eV (Sakurai et al. 2012). Strongly magnetized NSs, in- 
cluding Her X— 1 and GX 1+4, in contrast show narrow iron lines with significantly larger EWs 
(> 50 eV; Torrejon et al. 2010). These differences can be employed as an additional empirical 
argument to strengthen the high B scenario of 4U 1822—37, although theoretical account of 
this issue is beyond the scope of the present paper. In this respect, The narrow (cr=0.04 keV) 
and strong (EW~ 50 eV) Fe Ka line, detected at 6.4 keV with the XIS, clearly classifies this 
NS into the high-field category. 

To be somewhat more quantitative, the width of the Fe K a line, 0.04+0.02 keV, implies 
that the velocity v gas of the gases which emit the iron line is ~ 0.6% of the light speed. Assuming 
that the gas obeys Keplar rotation, its distance becomes r gas = GM /v gas ~5.7 x 10 9 cm {G being 
the gravitational constant and M ~ 1.4M 0 the NS mass), which is comparable to the Alfven 
radius (~ 10 8 ~ 9 cm) for B ~ 10 12 G while much larger than the NS radius (~ 10 6 cm). Therefore, 
the matter creating the fluorescent Fe-K line is likely to be stored on the Alfven surface. 

In the present data, we detected not only the nearly neutral Fe K a line, but also the 
other two lines which are identified as Helium like (Fe XXV) line and Hydrogen like iron (Fe 
XXVI) line as judged from their line energies. Empirically, such ionized iron lines are observed 
from X-ray pulsars mainly when the source luminosity is high; e.g. 5 Lx >lxl0 3 ' ergs 1 . Such 
examples include Cen X-3 (Naik et al. 2011), Her X-l (Ji et al. 2009), LMC X-4 (Neilsen 
et al. 2009), and Be pulsars in luminous outbursts (Naik et al. 2013). Quantitatively this is 
reasonable, because ionization of circum-source materials is determined by so-called ionization 
parameter £=Lx/{n e ■ r 2 ), where n e is the electron density of the line emitter while r is its 
distance from the NS (ionizing photon source). Assuming that the two ionized lines in 4U 
1822—37 are emitted by the same material that is different from those emitting the 6.4 keV 
line, we estimate as £ ~2000, from Kallman & McCray (1982) and the comparable EWs of the 
two ionized lines (tablel). Employing the same argument as for r gas , we may obtain r~3x 10 9 
cm from the line width of ~ 0.05 keV. Adopting n e = 5 x 10 15 cnr 3 from Iaria et al. (2013), the 
luminosity is lxlO 3 ' erg s _1 , in agreement with estimate using equation (2). Thus, the source 
is considered to be relatively luminous in X-rays. 

During the net exposure of 37 ks (~ lOh), no X-ray bursts were detected. Similarly, 
there have been not previous reports of burst detection from this objects. While this could be 
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due to rather infrequent burst occurrence at Lx > 2 x 10 37 erg s 1 (Cornelisse et al. 2003), a 
more appropriate explanation would be to regard 4U 1822—37 as an X-ray pulsar, which do 
not produce X-ray bursts. 

We hence conclude that 4U 1822—37 is another example of LMXB that contains a 
strongly magnetized ( B ~3xl0 12 G) NS. Furthermore, its luminosity is likely to be ~3xl0 3 ' 
erg s _1 , instead of the previously reported value of ~lxl0 36 erg s _1 , although it is at present 
unclear whether this discrepancy is due to inaccurate distance estimate, on due to the edge-on 
source geometry which could reduce the X-ray flux reaching us. 

We thank all members of the Suzaku hardware and software teams and the Science 
Working Group. M.S., K.M. and T.E. are supported by the Japan Society for the Promotion 
of Science (JSPS) Research Fellowship for Young Scientists, the Grant-in-Aid for Scientific 
Research (A) (23244024) from JSPS, and Grant-in-Aid for JSPS Fellows, 24-3320, respectively. 
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Table 1. Results of the model fit to the 1-50 keV Suzaku spectra of 4U 1822—37. 
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